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We introduce an integrated algorithm for full scale investigation of the earth regeneration effect 
of solar neutrinos with all three active neutrinos. We illustrate that the earth effect on low energy 
solar neutrinos (E v < 20 MeV) causes large variations in the survival probabilities for large values 
of the mixing angle #12 for a given value of #13. But it is maximal for small values of #13 and 
diminishes for large values at a given value of 812. The nadir angle dependence is sensitive to the 
value of E' = B/A12. As far as low energy solar neutrinos are concerned we notice that the earth 
effect is more pronounced for core penetration. We confirm that the earth effect leaves the SMA 
region almost intact regardless of #13 and that it is significant for high energy 8 B and hep neutrinos 
in the LMA region. We also show that for large values of 813, the earth effect is significant for high 
energy neutrinos (E' > 10 9 ) and it can be large at any nadir angle. 
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I. INTRODUCTION 



Neutrinos are produced abundantly from thermonuclear reactions taking place in the stellar core. They carry away 
I most of the energy during the supernova explosion. Due to the extreme weakness of its interaction it escapes its 
■ cradle almost intact and provides us with most reliable clues to the internal physical environments of a star. Thus 
neutrino astronomy has been established as an important research field recently. The Standard Solar model (SSM) 
jl]-[3| predicts the flux of the electron neutrino. The historic Homestake experiment B for detecting solar neutrinos 
revealed that the detected neutrinos fall far short of the SSM prediction. 

Several mechanisms were proposed for the depletion of solar neutrinos. Pontecorvo [EJ proposed the idea of neutrino 
flavor oscillation earlier than any unification theories (6). (see rcf. [0 for an interesting review). The possibility of 
I solar neutrino flux reduction by vacuum oscillation was investigated earlier |^||. Subsequent experiments jl^ |l3| ] 
^"i' revealed that the solar neutrino deficit depends on the energy of the neutrino. Recently the idea of flavor oscillation 



O !* combined with the Mikheyev-Smirnov-Wolfenstein (MSW) effect 14, T^] has been most popular. The MSW equations 



X 



have been vigorously studied. For simple matter densities, some authors |Tq-|20|| derived accurate analytic solutions 
' for the two-neutrino MSW evolution equations. These solutions have been used for computing the neutrino survival 
probabilities with the approximation of the solar matter density by a mosaic of exponential functions. Direct numerical 
analysis was carried out also [^l). Two- flavor MSW solutions to the solar neutrino puzzle were fully analysed with 
high accuracy, including both the day /ni ght effect and spectral information ]2^f23| ]. 

Recent Super-Kamiokande experiment [p4|| shows an evidence for oscillations of the atmospheric neutrinos [^5] Q . 
The data are in good agreement with two-flavor «-> v T oscillations. These results did establish that neutrinos 
oscillate and possess non-zero masses. Thus full three-flavor neutrino oscillation analysis is needed for more accurate 
estimation of neutrino masses and mixing anlges. Even four neutrino solutions with the sterile neutrino have been 
suggested Hf-fl- 

It has been argued |53] -Q that under certain conditions the solution of the solar neutrino evolution equation 
in the two neutrino case can be extended to the three neutrino case with a slight modification. There have been 
many analytic works on three-flavor MSW solutions of the solar neutrino problem with simple linear and exponential 
matter densities [p3|-40|. Most three-flavor MSW fits to the solar neutrino experiments have been made within the 
above framework |41 — |45f| . Recently exact solutions to the three neutrino MSW equations for such simple matter 
densities were derived 46-4^]. Direct numerical analysis is always desirable but it has been hampered by the huge 
computational costs. Recently we presented an efficient numerical algorithm for direct computation of the solar 
neutrino survival probability with all three active neutrinos BSpfJ. 
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According to the PREM model ^J, the earth consists of 10 layers of different matter densities and compositions. 
The PREM model gives the matter density in each earth layer as a polynomial in radius. It is around p — 1.6 ~ 6.1 
mol/cm 3 . In such matter densities, the MSW effect can be significant for neutrinos with certain energies. The earth 
regeneration effect for solar neutrinos was numerically analysed by many people within the two-flavor framework 
p2|-|63t. Since numerical analysis takes a lot of computing time a few attempts to derive analytic formula were made. 
Approximating the matter density in each of the five simplified earth layers as an even degree quartic polynonimal in 
radius, Lisi and Montanino derived analytic solutions for the two-neutrino MSW equations. Other attempts were 
made for three neutrino MSW equations with a periodic step-function density profile [p5| and with constant matter 



densities |66 . However, to the best of our knowledge, there are no direct numerical works on the earth regeneration 
effect with all three neutrinos to the depth as done in the two flavor case. 

In this paper we present an integrated algorithm for analyzing the earth regeneration effect of solar neutrinos 
with all three active neutrinos considered. Integration of neutrino evolution equations can be done with any ordinary 
differential equation solver such as adaptive Runge-Kutta or Bulirsch-Stoer methods (see, e.g. p7|). The wavefunction 
of a neutrino passing through the earth oscillates rapidly with large amplitudes in certain ranges of energy and mixing 
parameters. The crux of the problem is how to take time averages out of these wildly oscillating wavefunctions. Two 
methods have been used in the two neutrino case. Baltz and Weneser |54|] carefully analyzed the oscillatory behavior 
of the survival probabilities and derived a formula that uses instantaneous values of the two particular wavefunctions. 
Exploiting the symmetry of the transition probabilities, Smirnov [ |59| derived a formula that uses only one transition 
probability, P{yi — * v e ). Due to economy of the Smirnov's formula it has been widely used. However, we could not 
extend the formula to the three neutrino case. There is no shortcut in this case and we need to integrate the evolution 
equation for three different initial conditions, we realized. We introduce a method to obtain time averages of survival 
probabilities of solar neutrinos after passing through the earth in the three neutrino case. Our method is numerical 
and thus quite general and can be used with any number of flavors. 

This paper is organized as follows: In section 2 we review the basic formalism very briefly. In section 3 we introduce 
our method for taking averages of neutrino fluxes as they pass through the earth core. In section 4 we illustrate the 
complexity of the earth effect occurring at large mixing angles. We illustrate how it depends on neutrino energies, 
the mixing parameters, and the nadir angles. In section 5 we present how to make long term averages such as daily 
averages in different seasons and monthly and yearly averages. In section 6 we present sample contour plots for an 
interesting set of parameters. 



II. BASIC FORMALISM 



The state of a neutrino can be expressed either in the weak eigenstate basis \v a ) or in the mass eigenstate basis 
\vi). The unitary matrix U transforms the mass eigenstates into the weak eigenstates, 
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Vg. = U a m, a = e,/j,,T. (1) 

The Particle Data Group J68| adopts the convention for the mixing matrix, 

(C1C3 S1C3 S3 \ 

— S1C2 — C1S3S2 C1C2 — S1S3S2 C3S2 , (2) 
S1S2 — C1S3C2 —C1S2 — S1S3C2 C3C2 J 

where Ci = cos 9i and Si = sin#j and the three mixing angles B\ = #12, O2 = O23, and #3 = #13 roughly measure 
mixing between mass eigenstates (1-2), (2-3), and (1-3) respectively. We have neglected the CP violating phase, which 
is irrelevant in this problem. 

It is convenient to express the state of a neutrino as a mixture of three active neutrinos in the weak eigenstate basis, 

$(i) = a e \u e ) + a M (t)K) + a T (t)\v T ). (3) 
The mass matrix is non-diagonal in the weak eigenstate basis. It is obtained from the diagonal mass matrix, 

Mij=J2 U ik m l U kj- (4) 

fc 

The wavefunctions a a (t) of a neutrino in a medium obey the evolution equations, 
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Oe(t) \ -i / Mn+A c (i) Mi 
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4 I 0^(4) I =^ I M 21 M 22 i\/ 2 , | j „„(/) | . (:,) 

a T (t) / ^ \ M 3 i M 32 



where A c = 2\f2GpEN e is added due to the interaction with the charged currents via the MSW mechanism |i~4 15 
and (3 = E/R E (MeV/cm) and t = R/R E . 

In principle we can use the neutrino evolution equations in either of the two basis. In practice the mass eigenstate 
based equations require diagonalization of the mass matrix at each time step. Diagonalization introduces non-trivial 
errors and takes extra computing time in addition to the time to solve the ordinary differential equations. Thus we 
chose the weak eigenstate based evolution equations. 

Let us make a quick check of the two neutrino resonant condition in the earth core where N e ~ 0.466/9, 

cos(2%) = V2G F N e /(A i:j /2E), (6) 

while holding the other mixing angles equal to zero. One finds that for a small value of the (1-2) mixing angle 
(cos 20i2 ~ 1), the (1-2) transition takes place for neutrinos with E//S.12 ~ 10 6 ~ 10 7 eV _1 . For large values of 6*12, 
the resonant condition is satisfied for lower energy values. Thus the low energy pp-neutrinos can be susceptible to the 
earth matter. 



III. COMPUTATIONAL METHODS 



The electron neutrinos born from various sources within the solar core come out to the solar surface with average 
probabilities [J e 2 = P(v e — ► 24), X 2 = P(v e — > v^), X 2 = P(v e — > v T )]. In the interesting ranges of mixing 
parameters, A12, #12, 6*13, we have compiled a large set of solar neutrino survival probabilities born at selected radial 
positions r with selected energies E Q5(J. We used A 23 = 2.2 x 1CT 3 eV 2 and #23 = 43.5°. 

Integration of the neutrino evolution equations in the earth are carried out using the CERN library (69) routine 
DDEQBS, which is an implementation of the Bulirsch-Stoer algorithm. We used the PREM model fill with N e = 
0.466,0 in the core (R < 5463km) and N e = 0. 494,0 in the mantle. Our results agree with previous works |53] , p4] , |66[ . We 
especially checked our results against those of |36j and found good agreements, (see Fig. la). It would be extremely 
time-consuming to trace all these neutrinos individually as they pass through the earth. Instead we trace how three 
pure neutrinos, ^ e (t = 0) = (1,0,0,0,0,0), = 0) = (0,0,1,0,0,0) and ^ T {t = 0) = (0,0,0,0,1,0) evolve as 

they pass through the earth for a given set of parameters, Ai 2 , A 23 , #12, #13, #23- We compile the wavefunctions for 
selected values of the energy E and the zenith angle rj. Then we set up a handy interpolation function to be used for 
computing the neutrino fluxes arriving at various experiments in any given period of time. 

Let us denote the states of pure neutrinos after passing through the earth as, 

* e (t) = a e (t)\v e ) + a^(t)\^) + a T (t)K), 

%(t)=M*)\^)+0M v ^)+Mt)Wr), (7) 

*r(*) = 7eWK)+7/ 1 WI^)+7r(t)K)- 

Then the state of a solar neutrino in a mixed state with probabilities, (X^,X^, X^) can be written as a linear sum of 
three pure states, 

VseQ) = X e exp(i0 1 )* e (i) + X„ exp(^ 2 )* AI (t) + X T * T (t), (8) 

where we assumed some phases for solar neutrinos while keeping the phase of the r neutrino zero. 
Let us denote the state of a solar neutrino arriving at the earth as, 

# s (t = 0) = a e \v e ) + a^v^) + a T \v T ). (9) 

Then the probability of a solar neutrino to emerge as an electron neutrino after passing though the earth is written 

as 

Pse = |(^ e |*SB)| 2 = \a e u e + a-fiPe + a T ~f e \ 2 
= |a e | 2 K| 2 + | aAI | 2 |/3 e | 2 + K| 2 | 7e | 2 

+ [a e a e a*/3* + ex.] + [a^/3 e a*7* + c.c] + [a T j e a* e a* e + c.c.]. (10) 

In order to obtain the time average of PsE{t) it is more convenient to use the mass eigenstate basis. Let us denote 
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Oe(t) = sie 1 " 1 * + S2 e^ 2t + s 3 e l ^\ (11) 

/3 e (t) = t ie iuJlt + t 2 e luJ2t + t 3 e tuJ3t , (12) 

7e(t) = uie lWlt + u 2 e lW2t + w 3 e M3 ', (13) 

where Sj, tj, it, are independent of time. Then we have 

a e a* = Sis* + s 2 s 2 + S3S3 + cross terms, (14) 



a e [3* — s%tl + s 2 t 2 + 53*3 + cross terms 

= (sigtiR + s u t u ) + i{s u t 1R - si R t u ) + (s 2R t 2R + s 2I t 2I ) + i{s 2 it 2R - s 2R t 2I ) 

+ (s 3R t 3R + s 3I t 3I ) + i(s 3I t 3R - s 3R t 3I ) + cross terms (15) 

= Ra/3 + Haf3 + cross terms, 



where 



Rap = (smtiR + sut u ) + {s 2R t 2 R + s 2I t 2I ) + (s 3R t 3R + 837*37), (16) 
lap = (sutiR - si R tu) + (s 2 it 2R - s 2R t 2I ) + (s 3I t 3R - s 3R t 3 i). (17) 

The cross terms cancel out upon time averaging. 

It is convenient to express the amplitudes (a e , a^, a T ) in the weak eigenstate basis in terms of those in the mass 
eigenstate basis, (01,02,03), 

a e (t) = U n ai(t) + U 12 a 2 (t) + U 13 a 3 (t), 

a„(t) = U 2l ax{t) + U 22 a 2 (t) + U 23 a 3 (t), (18) 
Or(t) = U 3iai (t) + U 32 a 2 (t) + U 33 a 3 {t). 

Averaging over time, we obtain 

Xl = (a e a* e ) = UlMi? + Uf 2 \a 2 \ 2 + U 2 13 \a 3 \\- ■ ■. (19) 
We solve this equation for |ai| 2 . Then we can express (a e a*) in terms of | | 2 . Let us denote 



(a P a, 



UnU 2 iaial + U\ 2 U 22 a 2 a\ + U 13 U 23 a 3 a* 3 = X ell + iY efl . (20) 



We notice that the phases of become irrelevant upon time-averaging. Now we get the time average 

([a e a e a*/3* + c.c.]) = [(R a p + U a p){X ell + iY eiM ) + c.c] = 2{R a pX efl - I a pY ell ). (21) 

Similarly working on other terms and collecting all terms, we obtain the time averaged probability of a solar neutrino 
to emerge as an electron neutrino after passing the earth, 

(p S e) = kl 2 KI 2 + K\ 2 \p e \ 2 + kl 2 l7e| 2 

+ Z(R a pX etl — IapYefi + RfiyX^T- — Ip-yY^ T + R ia X Te — Iy a Y Te ). (22) 

In order to check our algorithm we reproduced some figures published in the literatures |54||34j in the limit 8\ 3 — > 0. 
(see Fig. 1). Our figures agree with the results as illustrated in Fig. 4 of j54j in the range E < 10 s . Our figures 
completely agree with the results obtained by using the time averaged Pe{v 2 — > v e ) in Eqn. (1) of [p4[. 



IV. SURVIVAL PROBABILITIES OF A SOLAR NEUTRINO WITH NON-ZERO 13 

We consider the effects of non-zero #13 on the electron neutrino survival probabilities of an electron neutrino 
created at the center of the Sun as a function of E' = E/&12 in MeV/eV 2 . We illustrate two cases in Figs. 2. In 
Fig. 2a we consider the Small Mixing Angle (SMA) case using the set of parameter values, (sin (2$i2) = 4.40 x 10~ 3 , 
A12 = 6.31 x 10~ 6 eV 2 ). It consists of 3 plots corresponding to #13 = 5°, 22.5°, 45°. Fig. 2b corresponds to the 
Large Mixing Angle (LMA) case and (sin 2 (2(?i2) = 0.76, A12 = 1.8 x 10~ 5 eV 2 ) was used. The upper limit of solar 
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neutrino energy, 20 MeV is indicated by the downward arrow. Each plot contains a curve for raw solar neutrino flux 
and five curves marked by Sm corresponding to five neutrino paths with nadir angles r\ — (m/10)(7r/2). 

Neutrinos undergo resonance conspicuously in those regions where the curves oscillate widely away from the raw 
curve. We notice that for nonzero values of the (1-3) mixing angle #13 double resonances occur, one in the range 
10 6 < E' < 10 7 and the other in 10 s < E' < 10 10 . The resonance in the low range is caused by non-zero #12 and its 
amplitude grows larger for large 6*12. The resonance in the high range is due to non-zero #13 and it is enhanced for 
large #13. The resonances occur for all nadir angles with almost equal strengths. It is noteworthy that for low energy 
solar neutrinos the earth effect is larger for smaller values of the nadir angle. 

Fig. 2a reveals that in the SMA case the earth effect for low energy solar neutrinos (E < 20 MeV) is minor for all 
values of 6*13. Fig. 2b exhibits that in the LMA case the earth effect is non-trivial for some high energy 8 B and hep 
solar neutrinos. It is most significant for small values of 6*13 and diminishes for large 6*13. 

Since it costs a lot of computation to trace neutrinos for each nadir angle we need to make a strategy for handling 
the angular effect. We exploit that for a given set of parameters, the neutrino wavefunctions are smooth functions of 
the nadir angle rj. We thus compute at a finite number of angles and then build an interpolation function. For the 
nadir angle, we choose cos(?7) as the sampling function. Thus we sample more frequently near rj = and scarcely near 
r\ = n/ 2. We selected 25 nadir angles between and tt/2. For low energy solar neutrinos the earth effect diminishes 
for 77 > r) c — 0.5779, the nadir angle grazing the boundary between the core and the mantle. Thus our sampling 
strategy works fine. An interpolation function made from these 25 rj values will be fairly accurate to use. 

We can now compute the states of pure neutrinos, ^ e {t), ^u(i), ^ T (t) at discrete energies and nadir angles selected 
according to the importance sampling functions. In our previous works on solar neutrino flux p9[ , we selected energy 
sampling points according to the convolution of the neutrino energy distribution functions as they are created in the 
solar core and the detector cross sections. In the present work we added 7 more energy sample points in the interval 
0.23 < E < 0.737 MeV to account for the CNO neutrinos more faithfully. 



V. LONG TERM AVERAGES: DAILY, MONTHLY AND YEARLY AVERAGES 

We have compiled a large amount of raw solar neutrino fluxes for a reasonably large area in the parameter space, 
(A12, sin 2 (26*12)) for selected values of #13, using the method described in p9] |. Our results were partially reported 
without the earth effect in [Q. The results of our computations in the 6*13 — > limit have been carefully checked 
against the results obtained with the analytic approximation formula such as the Petcov's improved formula Jl8[ . We 
agree with their results in the range of creation positions where their results are reliable as they claimed, in the sense 
that their formula give the flux values averaged over the creation positions whereas ours give the fluxes created at 
exact positions as illustrated in Fig. 5 of Ref. [Q. We also compared the survival probabilities of a neutrino created 
at the center of the Sun in the range of E 1 for various sets of mixing parameters. The agreement is to one part in 1000 
in the worst case. Since our results are computed numerically we have equally accurate fluxes beyond their claimed 
range. For three neutrino cases, analytic formula for survival probabilities for simple densities have been derived only 



recently [46-48J]. We have already confirmed that our subroutines reproduce the same results as those given in 

Next we compute seasonal variations of daily averages. The nadir angle rj of a neutrino path can be written |6j 
as a function of the date of year from the winter solstice, the time of day from the midnight and the lattitude 
A as follows, 

cos(?7) = cos(A) cos(r/ l ) cos(Js) — sin(A) sin(£s) (23) 

where the Sun declination is given by sin(i5s) = — sin(i£;) cos^) and %e = 0.4091 is the earth inclination. 

We illustrate the survival probabilities for two interesting parameter sets in Figs. 3. The average was taken over 
every 12 minuites for 24 hours on three days, the winter solstice (dotted), the vernal equinox (dash-dotted) and the 
summer solstice (dashed). The raw solar neutrino flux is plotted with the solid curve. The lattitude of Gran Sasso 
was used. The downward arrows mark the point of E = 20 MeV, the upper bound of the solar neutrino energy. Each 
set consists of 4 plots corresponding to #13 = 0°, 1.5°, 15°, 22.5°. As we notice from Fig. 3a, in the small mixing 
angle case the electron neutrino survival probability is almost unaffected by the earth effect as in the two neutrino 
case, regardless of 6*13. Fig. 3b shows that in the large mixing angle case (we used sin 2 (26'i2) = 0.63, A12 = 1.3 x 10~ 5 
eV 2 in Fig. 3b) the earth effect is non-negligible for high energy neutrinos. The effect is maximal for #13 = and 
diminishes as 6*13 is increased. 

We are now equipped with all the routines that are needed to compute average fluxes in an extended period of 
time. Since we used the importance sampling method for selecting reasonably large number of creation positions and 
energies to set up an interpolation function that is particularly good in the important ranges, we can reproduce the 
fluxes with high accuracy for creation position and energy values within the those ranges. 
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We have computed a large number of states of pure neutrinos, 9 e (t), ^r(t) at 100 discrete energies and 

25 nadir angles. We used 40 energy sample points for pp neutrinos in the interval 0.23 < E < 0.737 MeV similarly 
picked as in |^9|,^0| and the same 60 points as used in |49j5^| in the interval 0.737 < E < 19 MeV for boron 
and CNO neutrinos, for the same set of parameters (A12, sirr(20i2), 6*13). We handpicked three important energy 
points, E = 0.384, 0.862, 1.442 MeV for Be neutrinos and pep neutrinos repectively. Each integration of the 
neutrino evolution equations in the earth takes only a modest amount of computing time but the huge number of 
integrations requires a large amount of total computing time. We again have to resort to a parallel supercomputer. 
The parallelization can be done using the same method as described in |49|j . 

Using these data we can compute the earth regeneration effect on the survival probability of an incident solar 
neutrino at any energy and any nadir angle for a selected set of mixing parameters. In this way we can compute the 
earth regenerated solar neutrino fluxes for a large region of the mixing parameter space with a modest number of 
energy samplings to yield an accurate estimate of event rates. 

Some cautions are needed for setting up an interpolation function for creation positions and energies. We illustrate 
the solar neutrino survival probability as a function of creation position in Fig. 4, where we see a very sharply rising 
curve. We see a sharply dropping curve in Fig. 2a. Since we are using only a finite number of samples with different 
spacings we will encounter a sticky situation when the curves happen to turn sharply in the region where the sample 
spacing is wide. The popular cubic spline interpolation fails badly in this situation. One needs to use the exponential 
spline 0. We used the FITPACK subroutines available in the GAMS library Q. As Fig. 4 shows the interpolated 
curve is as good as the original one. 



VI. SAMPLE CONTOUR PLOTS 



We illustrate sample iso-SNU(FLUX) contour plots in the plane, (sin (2#i2) — A12), for a set of mixing parameters, 
6»i3 = 10°, A 23 = 2.2 x 10~ 3 eV 2 , and 6 23 = 43.5° in Fig. 5. For a given set of (sin 2 (26»i2), A i2 ) we had generated 
solar neutrinos according to the SSM model jl],H|. In order to compute the raw solar neutrino fluxes we used 36 
creation positions and 60 energy samples for pp neutrino and 60 creation positions and 60+7 energy samples for 
others. For the same set of mixing parameters (sin (2#i2), A12) we had computed the states of three pure neutrinos, 
^ e (t), *&fi(t)) ^t(U at 100 discrete energies and 25 nadir angles optimally selected as explained in section IV. For each 
solar neutrino created at position r with energy E we perform a sequence of time averaging operations leading to Eq. 
(22). In this way we build an interpolation function for (Pse{t, E, 77)) as a function of creation position r, energy E 
and nadir angle rj. We integrate the fluxes over the creation positions according to the relevant creation probability 
distribution functions and then build an interpolation function for (Pse{E, T))) as a function of the energy and nadir 
angle. We can now take seasonal and yearly averages of this flux to get (Pse(E)), which involves integration over i] 
with a definite lattitude A. Finally we take a convolution of (Pse{E)) with the absorption cross sections <j{E) to get 
the yearly average of event rate for a particular experiment. 

In order to make iso-SNU(FLUX) contour plots, we used a 16 x 21 logarithmic mesh in the (sin 2 (2^i2) — A12) 
plane in the range 0.001 < sin 2 (26>i 2 ) < 1 and 1.58 x 10 -6 < A 12 < 1.58 x 10~ 4 eV 2 . We padded two extra values, 
(0.85,0.97), for sin 2 (2#i2) in order to examine the large 612 region more carefully. 

We made separate plots for each experiment [p~0|— p~3|| with appropriate lattitude A. As expected from previous figures, 
contours in the SMA region are almost intact but those in the small A12 < 10~ 5 eV 2 and large sin 2 (2#i2) > 0.1 region 
are distorted significantly by the earth effect. But for gallium detectors the curves are not distorted visibly in the 
considered region of parameter space. This is due to the fact that the major signal source for gallium detectors is the 
pp neutrinos with very low energies, for which the earth regeneration effect is minimal. The overall earth effect for 
solar neutrinos is larger for small #13 than for large #13 as we see from Fig. 6. 



VII. SUMMARY AND CONCLUSION 



We introduced an integrated algorithm to deal with the earth regeneration effect of solar neutrinos with all three 
active neutrinos considered. Main ingredients of the algorithm are the time averaging algorithm with the use of mass 
eigenstates, the strategy for sampling energy and nadir angle, and the interpolation algorithm. Our algorithm is 
useful for full scale investigation of the earth effect. 

We illustrated that the earth effect on low energy solar neutrinos causes large variations in the survival probabilities 
for large mixing angles #12 at a given value of #13. But it is maximal for small values of #13 and diminishes for large 
values (~ 45°) at a given value of 6*12. The nadir angle dependence is sensitive to the value of E'. As far as low energy 
solar neutrinos are concerned we notice that the earth effect is more pronounced for smaller nadir angles. 
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We have shown that for large values of #13 , the earth effect is significant for high energy accelerator or atmospheric 
neutrinos (E' > 10 9 ) and it can be large at any nadir angle. 
This work was funded by POSTECH research fund. 
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Fig. 3a Daily averages of solar neutrino fluxes (SMA case) 

Solid curves are for raw flux, dotted ones for earth regenerated flux on the winter solstice, dashed 
ones on the summer solstice, dash-dotted ones on the equinox. All solar neutrinos are created 
at the center of the Sun. Downward arrows mark the 20 MeV upper limit of E v . We used the 
mixing parameters: Am| 3 = 2.2 x 10~ 3 eV 2 , #23 = 43.5°. The values of #13 are inscribed in the 
plots. For the SMA case: Am 2 2 = 6.31 x 10" 6 eV 2 , sin 2 (2#i 2 ) = 4.36 x 10" 3 . For the LMA case: 
Am 2 2 = 1.3 x 1CT 5 eV 2 , sin 2 (26»i 2 ) = 0.63. 
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Fig. 3b Daily averages of solar neutrino fluxes (LMA case) 

Solid curves are for raw flux, dotted ones for earth regenerated flux on the winter solstice, dashed 
ones on the summer solstice, dash-dotted ones on the equinox. All solar neutrinos are created 
at the center of the Sun. Downward arrows mark the 20 MeV upper limit of E v . We used the 
mixing parameters: Am| 3 = 2.2 x 10~ 3 eV 2 , #23 = 43.5°. The values of #13 are inscribed in the 
plots. For the SMA case: Am 2 2 = 6.31 x 10" 6 eV 2 , sin 2 (2#i 2 ) = 4.36 x 10" 3 . For the LMA case: 
Am 2 2 = 1.3 x 1CT 5 eV 2 , sin 2 (26»i 2 ) = 0.63. 
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Fig. 1 Comparison of earth regenerated solar neutrino fluxes: 

(top) Fig. (1) of Ohlsson Ref. [66], (middle) Fig. 3d of Baltz and Weneser (1987), Ref. [54] (bottom) Fig. 
5d of Ref. [54] 





Fig. 2a Earth regenerated solar neutrino fluxes (SMA case) 

The solid curve (DY) is for raw flux and the others marked by Sm are for the paths with the nadir angles 
7] = (to/10)(7t/2). All solar neutrinos are created at the center of the Sun. Downward arrows mark the 20 
MeV upper limit of E v . We used the mixing parameters: Am| 3 = 2.2 x 10~ 3 eV 2 , 023 = 43.5°. The values 
of 0i3 are inscribed in the plots. For the SMA case: Amf 2 = 6.31 x 1(T 6 eV 2 , sin 2 (20i 2 ) = 4.40 x 10~ 3 . 
For the LMA case: Amj 2 = 1.8 x 10~ B eV 2 , sin 2 (20i 2 ) = 0.76. 





Fig. 2b Earth regenerated solar neutrino fluxes (LMA case) 



The solid curve (DY) is for raw flux and the others marked by Sm are for the paths with the nadir angles 
7] = (to/10)(7t/2). All solar neutrinos are created at the center of the Sun. Downward arrows mark the 20 
MeV upper limit of E v . We used the mixing parameters: Am| 3 = 2.2 x 10~ 3 eV 2 , 023 = 43.5°. The values 
of 0i3 are inscribed in the plots. For the SMA case: Amf 2 = 6.31 x 1(T 6 eV 2 , sin 2 (20i 2 ) = 4.40 x 10~ 3 . 
For the LMA case: Amj 2 = 1.8 x 10~ B eV 2 , sin 2 (20i 2 ) = 0.76. 
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Fig. 4 Survival probability of a solar neutrino as a function of creation position 



Original data are marked by the solid curve, the cubic spline fit by the dashed and the exponential fit by 
the dotted. 




Fig. 5. iso-SNU(FLUX) plots for solar neutrinos with the CI, Ga, and Water detectors. 

The upper-left figure is for CL with SNU values 1.5, 2.56, 3.5, 4.5. The upper-right figure is for Water 
with FLUX values 1.5, 2.44, 3.5, 4.5. The lower-left figure is for GA with SNU values 55,65,73.1,90. 
Each contour is represented by dotted, dashed, long-dashed and solid curves respectively. In all these 
three figures the thick curves are drawn with the earth regenerated fluxes while the thin curves with raw 
fluxes. In each plot the vertical axis refers to Amf 2 in eV 2 and the horizontal axis to sin 2 (2#i2). We used 
Aml 3 = 2 - 2 x !0 -3 eV 2 , 8 23 = 43.5°, 8 13 = 10°. 




Fig. 6. iso-SNU(FLUX) plots for solar neutrinos with the CI detectors. 

Contours are for SNU values 1.5, 2.56, 3.5, 4.5. Each contour is represented by dotted, dashed, long- 
dashed and solid curves respectively. In all these four figures the thick curves are drawn with the earth 
regenerated fluxes while the thin curves with raw fluxes. In each plot the vertical axis refers to Afflj 2 in 
eV 2 and the horizontal axis to sin 2 (20i 2 ). We used Am| 3 = 2.2 x 1(T 3 eV 2 , 9 23 = 43.5°. 



